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ABSTRACT 

For a general mechanical system, it is shown that each solution of the 
Hamilton-Jacobi equation defines an N = 2 pseudo-supersymmetric ex- 
tension of the system, such that the usual relation of the momenta to 
Hamilton's principal function is the 'BPS' condition for preservation of 
1/2 pseudo-supersymmetry. The examples of the relativistic and non- 
relativistic particle, in a general potential, are worked through in detail 
and used to discuss the relation to cosmology and to supersymmetric quan- 
tum mechanics. 



1 Introduction 



The reparametrization invariant dynamics of a Hamiltonian system with s degrees of 
freedom is determined by a Hamiltonian constraint on a (2s + 2)-dimensional phase 
space. Given local Darboux coordinates (X m ; P m ) (m = 0, 1, . . . , s), the Lagrangian 
takes the form 

L = X m P m -£H(X;P) , (1.1) 

where the overdot indicates differentiation with respect to an arbitrary time parameter, 
and £ is a Lagrange multiplier for the Hamiltonian constraint. The Hamilton equations 
of motion are then 

p-lj C m_^L p-lp (12) 

The constraint function Ti will be assumed to be a polynomial in the momenta P that 
is at most quadratic. This is the case of most interest and most others can be put in 
this form by increasing the dimension of the phase space. Thus, 

U = X -K™\X)P m P n + J m (X)P m + U(X) , (1.3) 

for symmetric tensor field K{X), vector field J{X) and scalar potential function U(X). 
If K is invertible it may be interpreted as the metric on a target space with local 
coordinates X. The case in which K is not invertible 1 is also of interest because this 
allows for some components of P to appear linearly, but then we must insist that 
v ■ J 7^ for any co-vector field v(X) in the kernel of K; in other word 

v m J m ^0 if K mn v n = 0. (1.4) 

Otherwise, there are components of P that do not appear in Ti, so their conjugate 
variables are constants and we may rewrite the model in terms of a lower- dimensional 
phase space. 

The action on the constraint surface, viewed as a function of X at a 'final' time, is 
Hamilton's 'principal' function: 

S(X)= [ X dX m P m , (1.5) 



from which we we deduce that 

P m = d m S (m = 0,l,...,s), (1.6) 

where d m = d/dX m . These equations typically yield first-order differential equations 
for X when combined with the equation of motion for P, so we will sometimes refer 



lr The possibility that K(X) is invertible for some X and non-invertible for other X or, more 
generally, that the dimension of its kernel is position dependent, will not be considered here because 
neither it is considered in expositions of Hamilton- Jacobi theory. 
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to them as the 'first-order' equations of HJ theory. Using (1.6) in the Hamiltonian 
constraint, we find the Hamilton- Jacobi (HJ) equation 

^K mn d m Sd n S + J m d m S + U(X) = . (1.7) 

The solutions of this equation are in one-to-one correspondence with solutions of (1.2). 
This is not obvious, and Hamilton (who found the equation) did not appreciate this 
point, which is due to Jacobi. The main aim of this paper is to present a new deriva- 
tion of HJ theory, with various novel features, one of which is that it makes this 
feature of the HJ formalism manifest. The new derivation relies on a correspondence 
between solutions of the HJ equation for a given mechanical system and N = 2 pseudo- 
super symmetric extensions of it. 

The concept of pseudo-supersymmetry has its origins in supergravity [1, 2], and 
can be viewed as a complex analytic continuation of 'standard' supersymmetry. Its 
relevance to the Hamilton- Jacobi approach to inflationary cosmology was noticed in [3, 
4], where the 'first-order' equations of the HJ approach to cosmology [5] were recovered 
from an effective relativistic particle mechanics model, and interpreted as integrability 
conditions for the existence of 'pseudo-Killing' spinors 2 . This prompted several explicit 
realizations of pseudo-supersymmetric cosmologies in a ('variant') supergravity context, 
where the pseudo-Killing spinors become 'genuine' Killing spinors associated with the 
partial preservation of a symmetry [7, 8, 9]. In this sense, the equations (1.6) for these 
models may be interpreted as 'BPS' conditions, but this interpretation is not intrinsic 
to the effective particle mechanics model because it involves consideration of spacetime 
spinors. However, the supergravity examples do suggest the possibility of an intrinsic 
BPS interpretation, via an explicit pseudo-supersymmetric extension of the effective 
particle mechanics model. This was the starting point for the work reported here. 

It turns out that the required pseudo-supersymmetric extension of mechanics is 
not difficult to construct. As will be explained, there is a close relationship to stan- 
dard supersymmetric mechanics [10], but this suggests a difficulty: supersymmetry 
implies restrictions on the possible potentials, and one therefore expects the same to 
be true of pseudo-supersymmetry. In contrast, the HJ formalism does not involve any 
such restriction. Here, the key observation is that the superpotentials required for 
pseudo-supersymmetry may be multi- valued functions with branch points [3, 11]. In 
the context of the non-relativistic particle, the relation of this observation to HJ theory 
is obvious, and it leads directly to the general pseudo-supersymmetry formalism de- 
scribed here. We present this formalism from a classical perspective, initially, because 
the intention is to gain insight into standard results of classical mechanics, which we 
do by focusing on solutions of the equations of motion for which all anti-commuting 
variables are zero. 

In the quantum theory, one cannot set the anti-commuting variables to zero because 
this is not consistent with their anti-commutation relations. They could be 'integrated 

2 Because of the "Domain- Wall/Cosmology Correspondence" [6, 3], there is an analogous story for 
domain walls and many of the observations summarized here were first made in that context. However, 
in this introduction we focus on cosmology for simplicity of presentation. 
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out' (in a path-integral formulation) but this still implies some contribution to the 
wave-function. The standard semi-classical wavefunction has the form pexpiS for 
some variable modulus p determined by quantum fluctuations, so one would expect the 
anti-commuting variables to contribute to p. This is indeed the case, and the net result 
is that p is constant: the quantum fluctuations of the new anti-commuting variables 
cancel the quantum fluctuations of the original model. In effect, the quantum theory 
of the pseudo-supersymmetric mechanics model is equivalent to classical mechanics of 
the original model! This is reminiscent of the reformulation of classical mechanics of 
Gozzi et al. [12, 13] but the details appear to be rather different; in part because of 
the central role of the Hamilton- Jacobi equation in the formalism presented here. 

1.1 Particle Mechanics Examples 

As an aid to understanding the general formalism to follow, it may help to keep in 
mind the case of a particle in a space of dimension s, which yields various possible 
mechanical systems with s degrees of freedom. We shall later consider the s = 1 case 
for (i) a relativistic particle and (ii) a non-relativistic particle. When needed, we use 
the following notation for the components of the phase superspace coordinates: 

X m = (t,x), P m = (-E,p), 

A m = ftM), A m = HM). (1.8) 

The momentum P is a co- vector in a 2-dimensional Minkoski spacetime with signature 
(— 1, 1), and 

pi = r] mn P m P n = -E 1 + V 2 . (1.9) 
The two cases to be considered later are defined as follows: 

• Relativistic particle: 

h = \p 2 + u(x) . (i.io) 

This illustrates the case in which K is non-degenerate and J = 0. For 2U = m 2 we 
have a free relativistic particle of mass m, but we will consider a more complicated 
potential that is relevant to cosmology. 

• Non-relativistic particle: 

n = -E + ^ P 2 + u(t,x). (i.ii) 

This illustrates the case in which K is degenerate and J ^ 0. For t-independent 
potential we will write U — V(x). 
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2 The Formalism 



We begin by supposing that (X; P) phase space is the 'body' of a phase superspace 
with coordinates (X, A; P, A), where the anticommuting (s + l)-vector A and (s + 1)- 
covector A are canonically conjugate. We now allow for an extension of the Hamiltonian 
constraint function H on the (X; P) 'body' of the phase superspace to a function EI on 
the full phase superspace, and then consider a Lagrangian of the form 

L = X m P m + tA m A m -m + i X Q + ixQ (2-12) 

where (x, X) is a P arr °f re& l anticommuting Lagrange multipliers for a pair of con- 
straints with real anti-commuting constraint functions (Q, Q) . The Lagrangian itself is 
real because we adopt the convention that complex conjugation changes the order of 
anti-commuting quantities. It follows from this Lagrangian that the non-zero Poisson 
brackets of the dynamical variables are 

{X m ,P n }pB = — {P n , X m }pB = 8™ , 

{A m ,A n } PB = {A n , A m } PB = —i8™ , (2.13) 

and one may use this result to compute the Poisson brackets of the constraint functions, 
which must be in involution for consistency. We will require that the only non-zero 
Poisson bracket of constraint functions is 

{Q,Q}pb = -2iU. (2.14) 

If Q were the complex conjugate of Q then we would have a standard N = 2 super- 
symmetric extension of the model defined by (1.1). Instead, both Q and Q are real, so 
we have an N = 2 pseudo-supersymmetric extension. 

Each of the constraint functions generates a local symmetry. In particular, the 
infinitesimal pseudo-symmetry variations of any function <3> on the phase superspace 
are given by 

8 e $ = i {eQ, $} PB , = i {eQ, $} pB , (2.15) 

where the real anticommuting parameters (e, e) are arbitrary functions of the arbitrary 
independent variable. 

2.1 Principal function as Superpotential 

We shall now re-interpret the HJ equation (1.7) as an expression for U in terms of 
a superpotential S. We may then rewrite the Hamiltonian constraint function in the 
factorized form 

n = 1 [K mn + Q nS) + 2jm] _ dmS) (2 16) 

We now introduce a symmetric afline connexion V on the coordinate space such that 
both the tensor K and the vector J are covariantly constant: 

dpK mn = _ 2 r (m pq K n)q , d p J m = -T m pq J q . (2.17) 
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When K is invertible, and interpreted as a metric, the first of these conditions implies 
that T is the usual Levi-Civita connexion, and the second condition is then a constraint 
on the allowed choices for J. The curvature tensor constructed from T is 

R m npq = 2<9[ g r m p] n + 2T m e [ q T e p ] n , (2-18) 

and it satisfies the cyclic identity R m [ npq ] = 0. 
Now define 3 

F n = P n + tT m np A p A m , (2.19) 

and consider the choice 

Q = (F m -d m S)A m = (P m -d m S)A m 

Q = [K mn (F n + d n S) + 2J m )A m . (2.20) 

It is obvious that {Q, Q}pb = 0, and a computation of {Q, Q}pb shows that this is 
zero too provided that 

K p[r K nW R m] qep = 0, (2.21) 

which is a consequence of the cyclic identity for invertible K. Finally, a computation 
of {Q,Q}pb shows that (2.14) holds with 

H = ±[ K mn (F n + d n S) + 2J m ](F m -d m S) 

-iK m v (V p d n S) A n A m + 1 -K nl R m lpq A p A 9 A n A m . (2.22) 

As required, H — > H when all anti-commuting variables are set to zero. 

To summarize: for every solution S of the HJ equation of some given mechanical 
model, we have an N = 2 pseudo-supersymmetric extension of that model in which 
Hamilton's principal function S is re-interpreted as a superpotential. There is some- 
thing odd about this result: there was no restriction on the initial choice of potential 
U but we are now saying that it should be expressible in terms of a superpotential, 
so should this not restrict the potential in some way? The sharpest illustration of 
this 'paradox' is provided by the non-relativistic particle constrained to move on the 
x axis in a time-independent potential V(x). We shall be studying this example in 
detail; to anticipate, N = 2 pseudo-supersymmetry implies that V is given in terms of 
a superpotential W(x) by the formula 

V = E - l - (W) 2 , (2.23) 

where E is the particle's energy. Clearly, there are many potentials V that cannot be 
written in this form; for example, the harmonic oscillator potential V = x 2 . However, 
the potential V is only constrained if we assume that the superpotential W is defined for 

3 This shift of P 'covariantizes' the AA term in the Lagrangian, so all constraint functions become 
manifestly covariant when expressed in terms of P. 
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all values of x. There is no difficulty if we allow multi- valued superpotentials. In fact, 
the formula (2.23) is nothing other than the 'reduced' Hamilton- Jacobi equation for 
Hamilton's characteristic function W and, as is well known, the characteristic function 
has branch points at turning points of the motion. Thus, allowing for multi-valued 
superpotentials, a mechanical model has an N = 2 pseudo-supersymmetric extension 
for every solution of its HJ equation. 

2.2 The BPS condition 

Given a solution of the HJ equation, and hence an N = 2 pseudo-supersymmetric 
mechanics, we are faced with the problem of solving the constraints. One obvious way 
to do this is to set 

P = dS, A = 0. (2.24) 

Note that these equations are invariant under both pseudo-supersymmetries. This is 
obvious for Q and true for Q because of cancellations in the Poisson bracket with 
(P-dS). 

Another way to solve the constraints is to set 

K(P + dS) + 2J = 0, A = 0. (2.25) 

If K is non-invertible then this implies v m J m — for v in the kernel of K, which 
contradicts (1.4), so this alternative is viable only if K is invertible, and in this case it 
is equivalent to 

P + K^J + dS = 0, A = 0, (2.26) 

where 

P = P + K~ l J. (2.27) 
In terms of P the formula (2.16) becomes 

H=-(P + K~ X J + dS) k(P- K~ x J - dS) (2.28) 

from which we see that H. is invariant under P — > —P. This means that one gets 
equivalent physics by taking P — > —P, but making this transformation in (2.26) we 
recover the condition P = dS. 

Generically, there are no other ways to solve the constraints. For example, if one 
tries to consider a combination of the two alternatives just discussed, setting to zero 
some of the components of A and the complementary components of A, then one 
finds that the Hamiltonian constraint is not solved because there remain AA terms. 
This may not happen for special choices of S but then one expects additional symme- 
tries, which plausibly render any additional possibilities equivalent to (2.24). Given 
this, we conclude that the 'first-order' equations (1.6) are consequences of the pseudo- 
supersymmetry constraints. 
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Although, the conditions (2.24) preserve both pseudo-supersymmetries, a generic 
solution of the equations of motion for (X, A) will break both of them. To see this, 
first note that, 

5,A = 2 (KdS + J) e (2.29) 

for configurations satisfying (2.24), which shows that the Q pseudo-supersymmetry will 
be broken unless K8S+J = 0. This condition cannot be met when K is non-invertible, 
because of (1.4), so in this case the Q pseudo-supersymmetry is broken for all solutions 
of the equations of motion. Secondly, note that 

5 e X = -ieA , (2.30) 

which shows that the Q pseudo-supersymmetry will be broken unless A = 0. 

Our primary interest is in the original model, with (X; P) phase space, and solutions 
of this model are found from solutions of the pseudo-supersymmetric model by settting 

A = 0, A = 0. (2.31) 

Solutions of the equations of motion that have A = A = initially will have A = A = 
at all times, so we may consistently restrict attention to configurations of this type. In 
this case, the condition for partial preservation of pseudo-supersymmetry is that the 
variations of A and A vanish for non-zero e or non-zero e. The latter option is either 
not possible or equivalent to the former, so we assume that e = 0. The only non-zero 
pseudo-supersymmetry variation is then 

5 e A = (P-8S)e. (2.32) 

The condition for partial preservation of pseudo-supersymmetry is therefore P = 8S. 
The 'first-order' equation (1.6) of the HJ formalism may thus be interpreted as a 'BPS' 
condition for preservation of 1/2 pseudo-supersymmetry. Note, however, that this 
interpretation is simply a consequence of the constraints that define the model and the 
restriction to configurations satisfying (2.31): all solutions of the original mechanical 
model are BPS solutions of its pseudo-supersymmetric extension. 

2.3 Pseudo-Superspace 

Superfield methods may be used to make the pseudo-supersymmetries manifest. This is 
easily done for the Q pseudo-supersymmetry because the corresponding 'supercovariant 
derivative' D satisfies D 2 = 0, which can be realized as D = 8/89 for independent 
real anticommuting variable 9. We now interpret X and A as superfields with 9- 
components 4 

DX m = A m ? = • _ g mS ^ (2>33) 

4 As is customary, we use the same symbol to denote a superfield and its lowest component. No 
confusion arises as long as one employs a formalism (as we do) that allows any component equation 
to be viewed as a superfield equaton. 
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We also introduce the 'fermion number' superfield 



N = A m A. 



(2.34) 




(2.35) 



DQ 



2m. 



(2.36) 



Now consider the Lagrangian 



L = -iD X m L 



rn 



(2.37) 



The Lagrange multipliers £ and \ are now superfields with ^-components 



Dt = 2x, 



Dx = iq, 



(2.38) 



where q is a new commuting variable. One finds that, on omitting a total derivative, 
that 



where L is the Lagrangian of (2.12). The new variable q therefore imposes a constraint 
of vanishing 'fermion' number. This constraint is satisfied automatically by the solution 
(2.24) of the pseudo-supersymmetry constraints, which explains why we did not have 
to consider it previously. Also, it is consistent to set q = 0, in which case we recover 
(2.12) directly. 

We have now shown how the Q pseudo-supersymmetry may be made manifest. 
To do the same for the Q supersymmetry, the superfields X m and A m of the above 
discussion would have to be combined into a single N = 2 superfield depending on real 
anti-commuting variables 6 and 9, but this will not be attempted here, in part because 
it is not obvious how to proceed when K is non-invertible. 

2.4 Quantum Theory 

Let us now consider the quantum theory of the N = 2 pseudo-supersymmetric extension 
of some model of mechanics. In contrast to the viewpoint adopted so far, the additional 
anti-commuting variables cannot be 'set to zero' at the end because they are now 
operators acting on Hilbert space. In particular, this means that Q, Q and H become 
operators (denoted by a 'hat') that span a superalgebra for which the only non-zero 
(anti) commutator is 



L = L + iqN , 



(2.39) 




(2.40) 



As a consequence, we have the operator identities 



Q =0, 



(2.41) 
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For standard supersymmetry, Q is the hermitian adjoint of Q. For pseudo-supersymmetry, 
Q and Q are independent hermitian operators, given that EI is hermitian. However 
there are no nilpotent hermitian operators that act on a Hilbert space with positive 
definite norm. This might have been anticipated from the non-unitarity of pseudo- 
super symmetric field theories since one gets field equations from the first- quantization 
of particles. However, we set aside this difficulty for the moment. 

Upon quantization, the Poisson bracket relations of (2.14) become the (anti) com- 
mutation relations 



X m ,P n 



{A m ,A„} = C, 



which we may realize by 



X m = X m , A m = A m , P m = -id m , A m = d/dA r 
It may be verified that the following operators are nilpotent: 

Q = -iA m {d m -id rn S) , 



(2.42) 



(2.43) 



K mn (d n + id n S) + 2iJ m - T™Ki n J^A? 



d 



<9A r 



Nilpotency of Q is manifest and Q is nilpotent provided that 

K pq R lm piq K n]e = , 



(2.44) 



(2.45) 



which, for invertible K, is the statement that the Ricci tensor is symmetric. The 
operator H may now be defined by (2.40), and one finds that 



m = n + ... 

where the dots indicate terms that are annihilated by d/dA, and 

n = \ [K mn (p n + d n s) + 2 j m ] (p m - d m s 



(2.46) 



(2.47) 



The constraints are realized in the quantum theory by the following physical-state 
conditions on the wave-function ty(X,A): 



Q# = , = . 



(2.48) 



These are the only independent physical-state conditions since they imply that W^fr = 0. 
We solve them along the same lines as for the classical theory. Specifically, the classical 
conditions (2.24) become the following constraints on the wavefunction: 



d m V = i (d m S) V , d^/dA m = , 



for which the solution is 



m = e iS{x) ^ ( 



(2.49) 
(2.50) 
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for constant ^o- This is a rather surprising result because it implies that \I/ has constant 
modulus; this is possible because H is not hermitian, despite the (formal) hermiticity 
ofH. 

Another way of stating the above result is to observe that the effective action, after 
inclusion of quantum effects, equals the classical action, which means that the anti- 
commuting variables cancel out the quantum fluctuations of the original variables. In 
effect, they allow a quantum description of a classical theory. This is what is also 
achieved by the formalism of Gozzi et al. [12, 13] which involves additional anti- 
commuting variables that are interpreted as ghosts and anti-ghosts associated with a 
BRST and anti-BRST invariance. The phase superspace in that formalism has dimen- 
sion (4s 1 4s) rather than (2s|2s), and there is no special role for the HJ equation, so it is 
unclear whether there is any connection to the formalism presented here. Nevertheless, 
the possibility of such a connection suggests that the pseudo-supersymmetry charges 
introduced here should be interpreted as BRST charges. In this case the superspace 
expression (2.37) shows that the Lagrangian is BRST exact, and hence that we are 
dealing with a 'topological' theory, as argued in [14] for the formalism of [12, 13]. 

3 Relativistic Particle 

Let 

Q=(P m -d m S)A m , Q=(P m + d m S)A m , (3.1) 

where A m = i] mn A n . The only non-zero Poisson bracket of these functions is that of 
(2.14) provided that we choose 

H = l -rT n (P m + d m S) (P n - d n S) - i (d m d n S) A m A n . (3.2) 
This corresponds to a potential given by 

U = ~\ (dS) 2 , (3.3) 

which is just the HJ equation for our model. 

We may solve the pseudo-supersymmetry constraints in two equivalent ways: 

either P = OS , A = , or P = -OS , A = . (3.4) 

For either choice we may restrict to solutions of the equations of motion for which 
both A = and A = 0, and in this case either the Q or the Q pseudo-supersymmetry 
will be preserved, according to whether P = dS or P = —dS. The two possibilities 
are physically equivalent. Both pseudo-supersymmetries are preserved if (in addition) 
P = 0, but this is the 2-momentum of the vacuum. 
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3.1 Application to Cosmology 

The assumptions of homogeneity and isotropy in a model of gravity coupled to scalar 
fields with potential V, in D spacetime dimensions, lead to an effective relativistic 
particle mechanics model with a Hamiltonian constraint of the form (1.10). In the 
conventions of [3], the scale factor a is written as exp((3t), and the scalar potential of 
the effective particle mechanics model is 

U(t,x) = e 2at V(x)-^ a , (3.5) 

where k — — 1, 0, 1 is the normalized curvature of spatial sections, and 

a = (D-l)P, P=l/^2(D-l)(D-2). (3.6) 

There are two cases in which the resulting HJ equation may be solved by a separation 
of variables. The simplest is V = 0, in which case 

S = ± (D - 1) e t/{2a) [e x/{2a) - ke~ x/{2a) ] (V = 0). (3.7) 

The other case for which the HJ equation may be solved by separation of variables 
is k = 0, for which 

S = ±2e at W{x) = 2a D - 1 W{x) (k = 0) , (3.8) 

where the 'comoving principal function' W is t-independent and satisfies the 'reduced 
HJ equation' [5] 

V = -2 [(W'f - a 2 W 2 ] . [Cosmology] . (3.9) 
The equations (1.6) are now 

E = T2ae at W(x) , p = ±2e at W'(x) . (3.10) 

Combining these with the equations of motion for E and p we arrive at the following 
first-order differential equations: 

i~H = -2ae at W(x) , t x x = 2e at W\x) . (3.11) 

We have focused on cosmology but exactly the same analysis applies to domain 
walls. The effective action is the same, except for a flip of the sign of the potential, 
so an application of HJ theory to flat domain walls [15] leads again to the first-order 
equations (3.11) but with a 'reduced HJ equation' of opposite sign for V: 



V 



{W') 2 - a 2 W 2 . [Domain Wall] (3.12) 



Remarkably, this formula coincides with the 'supergravity-inspired' formula for V in- 
troduced in [16, 17]; in fact, for D = 3 it is the supergravity formula for V in terms 
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of a superpotential W. In the context of supergravity domain walls, first-order equa- 
tions consistent with (3.11) arise as 'BPS-type' conditions for the preservation of 1/2 
super symmetry [18, 19], and they can be found precisely in the HJ form (3.11) by 
'supergravity-inspired' methods [20, 21]. Specifically, they arise as integrability condi- 
tions for the existence of Killing spinors defined in the context of 'fake supergravity' 
[22, 23]. It was pointed out in [3] that this 'coincidence' between the HJ and fake su- 
pergravity approaches to domain walls also applies to cosmology but the (fake) Killing 
spinors become (fake) pseudo-Killing spinors, defined by an analytic continuation for 
which W — > iW . As mentioned in the introduction, it has been shown recently that 
fake pseudo-Killing spinors may be 'genuine' Killing spinors of a 'pseudo-supergravity' 
theory, defined by analytic continuation of a standard supergravity theory to one with 
spinors obeying 'twisted reality' conditions. 



3.1.1 Non-flat cosmology 

It was shown in [3] that non-flat cosmologies are determined by the following first-order 
equations involving a complex superpotential Z: 



E = ^2ae 



at 



Tie (ZZ') 



V 



±e at \Z'\ 



(3.13) 



These equation reduce to those of (3.10) when Z — W , which applies when k — 0. 
One may ask what the relation of these equations is to the first-order equations of HJ 
theory, 

E = -dtS, p = d x S. (3.14) 
If we try to combine (3.14) with (3.13), we arrive at the equation 



±dS = uj = 2ae 



at 



Tie ZZ' 



dt + e at \Z'\dx. 



(3.15) 



However, the 1-form uj is not closed for k ^ 0; using eq, (4.16) of [4], one finds that 



duj = — \Z'\- l e^ 2 ^dtNdx. 
2(3 1 1 



(3.16) 



It follows that (3.15) cannot be integrated to give a function S(t,x). An expression for 
S for k 7^ was found in [4] , but the construction assumed that one is given a solution 
of the equations of motion, in which case the variables (t,x) are not independent. 

For k — 1 cosmologies, the 'first-order equations that follow from (3.13) on using the 
equations of motion for E and p are integrability conditions for the existence of pseudo- 
Killing spinors, and in this 'field theoretic' sense they may be interpreted as 'BPS' 
conditions 5 , as for k — 0. Again, there is a parallel story for domain walls, most of which 
came first; we focus here on cosmology for convenience of presentation, but there is one 



'This property distinguishes these 'first-order equations from others proposed previously [24]. 
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point that is easier to understand from the domain-wall perspective: we do not expect 
to find Killing spinors for 'de Sitter sliced' walls, and this corresponds to the statement 
that we should not expect to find pseudo-Killing spinors for k = — 1 cosmologies. This 
is indeed the case, even though there are first order equations for any k This state 
of affairs should be contrasted with the particle mechanics BPS interpretation of the 
equations (3.14), which applies for all k. Moreover, in the particle mechanics BPS 
interpretation, the superpotential is real for all k, being identified with the 'comoving 
principal function' W. 

These comments suggest that there is no simple general connection between the field 
theoretic BPS interpretation of the first-order equations for cosmology that arise from 
the existence of pseudo-Killing spinors and the 'intrinsic' BPS interpretation proposed 
here in the context of an N = 2 pseudo-supersymmetric extension of mechanics. 



4 Non-relativistic particle 

For a non-relativistic particle in a potential U(t,x), the pseudo-supersymmetric exten- 
sion is found by choosing 6 

Q =( p -3 A+ ( s+ !h M p+ 3 x+2 *'- (4 - i7) 

where S(t,x) is a superpotential. The only non-zero Poisson bracket is that of (2.14) 
provided that we also choose 

B = -E + ^ + U-i(g)X + i(™y, (418, 

where 

which is the HJ equation of our model. 
To satisfy the constraints we set 

dS dS , inn . 

and 

^ = 0, A = 0. (4.21) 
From the infinitesimal pseudo-supersymmetry transformations 

Sip = -2e, 5X = (p+^j 1 



°This is suggested by a 'worldline supergravity' Lagrangian proposed in [25]. 
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we see that the Q pseudo-supersymmetry is necessarily broken. Provided that ip = 
A = 0, the Q pseudo-supersymmetry is preserved. 

4.1 Time-independent potential 

It will prove instructive to analyse in more detail the special case for which 

S(t,x) = W(x) - E t, (4.23) 

for constant E , and a function W that can be interpreted as Hamilton's "characteris- 
tic" function. In this case we have the t-independent Hamiltonian constraint function 

H = -E + \p 2 + V + iW"\\ , (4.24) 

where 

V = E - i (W'f , (4.25) 

which is the 'reduced' HJ equation appropriate for a principal function of the assumed 
form. The pseudo-supersymmetry constraint functions are now 

Q = (p - W') A + (E - E Q ) i> , Q = (p + W) A + 2^ , (4.26) 

and they generate the following infinitesimal pseudo-supersymmetry transformations: 

5x = -i (eA + eA) , 8p = -i (eA - eA) W" 

St = ieip, 5E = 

S\ = (p + w')e, 5X = (p-W')e 

dip = -2e, 5ip = -(E-E )e. (4.27) 

The Lagrangian (2.12) transforms into a total derivative under these transformations 
provided the Lagrange multipliers are assigned the transformations 

Sx = e, 5x = X, 5£=-2t(ex + e X ) ■ (4.28) 

To simplify things we will now partially fix the local pseudo-supersymmetry trans- 
formations by the gauge choice 

X = 0, x = 0. (4.29) 

This leaves a residual invariance under the transformations of (4.27) for constant (e, e). 
We must therefore still impose the constraints Q = and Q = at some initial time, 
but they will then hold at all times as a consequence of the equations of motion. The 
advantage of this gauge choice is that we now have the much simpler Lagrangian. 



L = xp + i\\ — £ 



E + ^p 2 -^(W') 2 + iW"\\ 



— iipip + E (i — tj . (4.30) 
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We may further simplify by using the equations of motion of the conjugate pair (if;, if;), 
which imply that 

if; = ip , ^ = -0o, (4-31) 

for anti-commuting constants (if) , if> ), together with the equations of motion of the 
conjugate pairs (t, —E): 

E = 0, £ = i. (4.32) 

Ultimately, one finds that the constant value of E must be Eq and we shall now assume 
this in order to shorten the presentation. We now have the even simpler Lagrangian 

L = xp + iXX-tW, (4.33) 

where, from (4.24) and (4.25), 

U = \p 2 -\ (W'f + iW'XX . (4.34) 

Z Z 

The pseudo-supersymmetry Noether charges are now 

Q = (p-W')X, Q= (p + W')X + 24> , (4.35) 



and they generate the transformations 7 



(4.36) 

(4.37) 
gauge choice 

i = 1 (4.38) 



5x = -i (eX + eAj , Sp = -i (eX - eAj W" , 
SX = (p + W')e, 5X = (p-W')e. 

together with 

5ip = -2e. 

Finally, we may fix the time-reparametrization invariance by the 



to arrive at the Lagrangian 

L = xp + iXX -H, (4.39) 

where H, which is still given by (4.34), may now be interpreted as the Hamiltonian. 
The gauge choice is again only partial because it leaves a residual invariance under rigid 
time translations, and this means that we must set H = at some initial time, now as 
a zero-charge condition which holds at all times as a consequence of the equations of 
motion. The zero charge conditions may be solved in the same way that we previously 
solved the phase superspace constraints. We must arrange for Q = and Q = in 
such a way that one of the two pseudo-supersymmetries is preserved when all anti- 
commuting variables are zero. 

7 The transformation St = ietpo is not generated by these simplified Noether charges because of the 
simplification we made in eliminating the variable E, but this is irrelevant because there is a manifest 
independent invariance under constant shifts of t that we may combine with the pseudo-supersymmetry 
to arrange for t to be inert. 



15 



Although we could omit the transformation (4.37) on the grounds that this is an 
independent symmetry (trivially since the Lagrangian is independent of tp ) we will 
retain it for the moment as a relic of the asymmetry between Q and Q that is inherent 
in the physics of the non-relativistic particle: this asymmetry is due to a choice of 
positive rather than negative energy in taking the non-relativistic limit. In this case, 
the symmetry generated by Q is necessarily broken, so we must preserve the symmetry 
generated by Q, and this requires 

p = W', A = . (4.40) 

We see that the equation relating the momentum to Hamilton's characteristic function 
is a BPS condition for 1/2 pseudo-supersymmetry. This should be combined with the 
equation (4.25) for V in terms of the superpotential W, which we may rewrite as 

W'(x) = ±yfe (E -V(x)) (4.41) 

and interpret as the reduced HJ equation for Hamilton's characteristic function W in 
terms of V. There is a solution only for Eq > V so if V is unbounded from above there 
is no choice of E for which W(x) is defined for all x. This is of course a well known 
result in mechanics. The corresponding result for superpotentials, and its implications 
in the context of supersymmetric domain walls was recently discussed in [11]. 

4.1.1 Relation to supersymmetric mechanics 

If we eliminate p from the Lagrangian (4.39) by using its equation of motion p — x, 
then we arrive at the equivalent Lagrangian 

L = -x 2 + iXX + - {W'f + iW'XX , (4.42) 
2 2 

which is invariant under the N = 2 pseudo-supersymmetry transformations 

Sx = -i (e\ + eX) , 5X = (x + W) e , 5X = (x-W')e. (4.43) 

This is similar to a model of N = 2 supersymmetric mechanics, but with 'wrong sign' 
potential, In fact, it is related to N = 2 supersymmetric mechanics by complex analytic 
continuation, as we now explain. 

Complexify all the dynamical variables (x,X, X). The result is a complex La- 
grangian, but one that depends analytically on them, and it remains invariant under 
the transformations (4.43) with complixified (e, e), which are now analytic transforma- 
tions. Clearly, we recover the original Lagrangian on restricting all fields to be real, but 
we may now ask whether there is some other way to get a real Lagrangian 8 . If there 
is, it must involve an analytic continuation of the real variables of (4.42). Consider the 
analytic continuation that effects 

W — > iW . (4.44) 
8 Here we follow the logic presented in [7] for supergravity. 
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This yields the new Lagrangian 

L = \x 2 + f\\-\(W'f - W"W. (4.45) 

This is again real provided we choose A to be complex with complex conjugate A. The 
same analytical continuation of the transformations (4.43) yields 

Sx = -i (eA + eA) , 5X = (x + iW) e , 5\ = (x- iW) e . (4.46) 

where e is now a complex parameter with complex conjugate e. This is now a standard 
N = 2 supersymmetric mechanics model with superpotential W. 

5 Discussion 

We have derived the Hamilton- Jacobi formulation of mechanics by considering the 
possible N = 2 locally pseudo-supersymmetric extensions of reparametrization invari- 
ant Lagrangians for mechanical systems for which the Hamiltonian is no more than 
quadratic in momenta, as can usually be arranged. The Hamilton- Jacobi equation 
arises as a condition for pseudo-supersymmetrizability, with the superpotential taking 
the role of Hamilton's principal function; one need not worry about the fact that there 
may not always be a solution because the superpotential is allowed to be multi- valued, 
with branch points, and there is no superpotential only when there is no real solution 
of the Hamilton- Jacobi equation. 

The 'first-order' equations of the HJ formalism relating the momenta to Hamil- 
ton's principal function arise as solutions to the constraints associated to the local 
symmetries, and if one sets all anti-commuting variables to zero then these constraints 
may be interpreted as BPS conditions arising from the preservation of 1/2 pseudo- 
super symmetry. This is very likely related to a 'field theoretic' BPS interpretation 
that arises from consideration of pseudo-Killing spinors in cosmological spacetimes for 
flat universes, which can be realized in some cases as supersymmetric solutions of 
pseudo-supergravity theories. However, a comparison with results of this 'field theo- 
retic' BPS interpretation for non-flat cosmologies makes it appear unlikely that there 
is any simple general relation. One possible route to a further investigation of this 
point would be to extend the considerations of this paper to mechanical models with 
an infinite number of degrees of freedom. After all, a field theory can be viewed as a 
model of mechanics on an infinite-dimesional space. From this perspective, the results 
obtained here should also apply to field theory, although functions such as Hamilton's 
principal function then become functionals. 

It was stated in the introduction that the pseudo-supersymmetric formulation of 
HJ theory makes it clear why the HJ equation is equivalent, if taken together with the 
'first-order' equations that define the momenta, to the Hamilton equations of motion. 
The reason is that in a reparametrization-invariant formalism, as used here, all the 
dynamics is encoded in the constraints. Thus, the dynamics of the model defined by 
(1.1) is encoded in the Hamiltonian constraint function H. However, to get to the HJ 
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equation we had to use the relations (1.6) for the momenta in terms of the principal 
function, and this needs a separate motivation (e.g. involving considerations of canon- 
ical transformations). In contrast, these relations are implied by the constraints in the 
pseudo-supersymmetrized model, while the model itself is determined by a solution 
of the HJ equations. Solutions of the original model are then found as solutions of 
the pseudo-supersymmetry constraints with vanishing anti-commuting variables, all of 
which preserve (at least) 1/2 of the pseudo-supersymmetry. We thus have a corre- 
spondence between solutions of the Hamilton equations (1.2) and solutions obtained 
by integration of the 'BPS' equations, alias the 'first-order' equations (1.6), for a given 
a solution of the HJ equation. 

Classical mechanics has been the focus of this work but it is natural to ask whether 
the anti-commuting variables arising from pseudo-supersymmetry have a role to play 
in the quantum theory. On the one hand, this seems to ruled out by the impossibility 
of a realization of the quantum pseudo-symmetry algebra via operators acting on a 
Hilbert space with positive definite norm. On the other hand, a simple computation of 
the wavefunction shows that the effect of the anti-commuting variables is to cancel the 
quantum fluctuations of the original variables, so the model is effectively still classical. 
This suggests a possible BRST interpretation of the pseudo-supersymmetry algebra, 
as in the formalism of Gozzi et al. [12, 13], which has some features in common with 
the formalism presented here. It would be satisfying if the two formalisms could be 
related; in any case, a better understanding of the quantum theory is clearly desirable. 
In view of the many insights into quantum mechanics provided by supersymmetric 
quantum mechanics (see e.g. [26]) it seems reasonable to hope that quantum pseudo- 
super symmetric mechanics will be useful too. 
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